On the stability of reduced-order linearized computational fluid dynamics models based on POD and Galerkin projection: descriptor vs non-descriptor forms

نویسندگان

  • David Amsallem
  • Charbel Farhat
چکیده

The Galerkin projection method based on modes generated by the Proper Orthogonal Decomposition (POD) technique is very popular for the dimensional reduction of linearized Computational Fluid Dynamics (CFD) models, among many other typically high-dimensional models in computational engineering. This, despite the fact that it cannot guarantee neither the optimality nor the stability of the Reduced-Order Models (ROMs) it constructs. Short of proposing any variant of this model order reduction method that guarantees the stability of its outcome, this paper contributes a best practice to its application to the construction of linearized CFD ROMs. It begins by establishing that whereas the solution snapshots computed using the descriptor and non-descriptor forms of the discretized Euler or Navier-Stokes equations are identical, the ROMs obtained by reducing these two alternative forms of the governing equations of interest are different. Focusing next on compressible fluid-structure interaction problems associated with computational aeroelasticity, this paper shows numerically that the POD-based fluid ROMs originating from the non-descriptor form of the governing linearized CFD equations tend to be unstable, but their counterparts originating from the descriptor form of these equations are typically stable and reliable for aeroelastic applications. Therefore, this paper argues that whereas many computations are performed in CFD codes using the non-descriptor form of discretized Euler and/or Navier-Stokes equations, POD-based model reduction in these codes should be performed using the descriptor form of these equations. David Amsallem Department of Aeronautics and Astronautics, Durand Building, 496 Lomita Mall, Stanford University, Stanford, 94305-4035, USA, e-mail: [email protected] Charbel Farhat Department of Aeronautics and Astronautics, Department of Mechanical Engineering, and Institute for Computational & Mathematical Engineering, Durand Building, 496 Lomita Mall, Stanford University, Stanford, 94305-4035, USA e-mail: [email protected]

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Extension Ability of Reduced Order Model of Unsteady Incompressible Flows Using a Combination of POD and Fourier Modes

In this article, an improved reduced order modelling approach, based on the proper orthogonal decomposition (POD) method, is presented. After projecting the governing equations of flow dynamics along the POD modes, a dynamical system was obtained. Normally, the classical reduced order models do not predict accurate time variations of flow variables due to some reasons. The response of the dynam...

متن کامل

Numerical Simulation and Parametric Reduced Order Modeling of the Natural Convection of Water-Copper Nanofluid

In this article, a coupled computational framework is presented for the numerical simulation of mass transfer under the effects of natural convection phenomena in a field contains water-copper Nano-fluid. This CFD model is build up based on accurate algorithms for spatial derivatives and time integration. The spatial derivatives have been calculated using first order upwind and second order cen...

متن کامل

A novel topological descriptor based on the expanded wiener index: Applications to QSPR/QSAR studies

In this paper, a novel topological index, named M-index, is introduced based on expanded form of the Wiener matrix. For constructing this index the atomic characteristics and the interaction of the vertices in a molecule are taken into account. The usefulness of the M-index is demonstrated by several QSPR/QSAR models for different physico-chemical properties and biological activities of a large...

متن کامل

Model Reduction of the Nonlinear Complex Ginzburg-Landau Equation

Reduced-order models of the nonlinear Complex Ginzburg-Landau (CGL) equation are computed using a nonlinear generalization of balanced truncation. The method involves Galerkin projection of the nonlinear dynamics onto modes determined by balanced truncation of a linearized system, and is compared to a standard method using projection onto Proper Orthogonal Decomposition (POD) modes computed fro...

متن کامل

Reduced order models based on local POD plus Galerkin projection

A method is presented to accelerate numerical simulations on parabolic problems using a numerical code and a Galerkin system (obtained vía POD plus Galerkin projection) on a sequence of interspersed intervals. The lengths of these intervals are chosen according to several basic ideas that include an a priori estímate of the error of the Galerkin approximation. Several improvements are introduce...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 2013